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ABSTRACT 

Emission lines of Be-like ions are frequently observed in astrophysical plasmas, and 
many are useful for density and temperature diagnostics. However, accurate atomic 
data for energy levels, radiative rates (A-values) and effective electron excitation colli¬ 
sion strengths (T) are required for reliable plasma modelling. In general it is reasonably 
straightforward to calculate energy levels and A- values to a high level of accuracy. By 
contrast, considerable effort is required to calculate T, and hence it is not always pos¬ 
sible to assess the accuracy of available data. Recently, two independent calculations 
(adopting the A-matrix method) but with different approaches (DARC and ICFT) 
have appeared for a range of Be-like ions. Therefore, in this work we compare the two 
sets of T, highlight the large discrepancies for a significant number of transitions and 
suggest possible reasons for these. 
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1 INTRODUCTION 

Emission lines of Be-like ions are widely detected in a va¬ 
riety of astrophysical plasmas, including the solar transi¬ 
tion region and corona, with the prominent ions including 
C III, O V, Ca XVII and Fe XXIII - see for example, the 
CHIANTI database at http://www.chiantidatabase.org/ 
and the Atomic Line List (v2.04) of Peter van Hoof at 
http://www.pa.uky.edu/~peter/atomic/. Many of the ob¬ 
served lines are sensitive to variations in density or tem¬ 
perature, and hence are useful as diagnostics - see for ex¬ 
ample, Landi et al. (2001). However, for the reliable mod¬ 


elling of plasmas, accurate atomic data are required, partic¬ 
ularly for energy levels, radiative rates (A-values), and exci¬ 
tation rates or equivalently the effective collision strengths 
(T), which are obtained from the electron impact collision 
strengths (H). Unfortunately, existing atomic data for Be¬ 
like ions prior to 2014 were very limited, particularly for T. 
However, very recently |Fernandez-Menchero, Del Zanna fcj 
Badnell (2014) have reported data for a range of Be-like ions 


up to Z = 36. For their calculations they have adopted the 
AutoStructure (AS) code of Badnell (1997) for the genera¬ 
tion of wave functions, i.e. to determine energy levels and 
A-values. In the subsequent calculations of H and T they 
have adopted the A-matrix code of |Berrington, Eissner fcj 
Norrington (19951. Their H are primarily obtained in the 


LS coupling (Russell-Saunders or spin-orbit coupling) and 
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corresponding results for fine-structure transitions are de¬ 
termined through their intermediate coupling frame trans¬ 
formation (ICFT) method. 

Assessing atomic data, particularly for T, is a difficult 
task (Aggarwal & Keenan 2013a I mainly because the calcu¬ 
lations are highly complicated, require large computational 
resources and cannot be easily repeated. Therefore, the as¬ 
sessment is frequently based on limited comparisons either 
with already available results or with simple calculations, 
such as with the distorted-wave method but without reso¬ 
nances. However, often such accuracy assessments are far 


from satisfactory as noted by Aggarwal & Keenan (2013a I 
for several ions. 

Realising the importance of Be-like ions we have also 


wal & Keenan 


reported atomic data for several, namely Al X (Aggar- 
2014a] ) , Cl X IV, K XVI a nd Ge XXIX 
& Keenanj| 2014c| ) and Ti XIX ( Aggarwal fe| 


Keenan 2012b I. Our calculations are independent of those 
of Fernandez-Menchero et ah] ([2014I and are based on the 
Dirac Atomic R-matrix Code (darc) of P. H. Norrington 
and I. P. Grant (http://web.am.qub.ac.uk/DARC/). Unlike 
the semi-relativistic version of the standard A-matrix code 


(Berrington et al. 1995) adopted by Fernandez-Menchero 


et al. (2014), DARC is based on the jj coupling scheme. 


The accuracy of the data calculated (for H and subsequently 
T) through this approach is generally higher, because reso¬ 
nances through the energies of degenerating levels are also 
taken into account. For this reason it particularly affects 
transitions among the fine-structure levels of a state. The 
degeneracy among such levels increases with increasing Z 
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Table 1 . Energies for the lowest 80 levels of A1 X (in Ryd). 


Index 

Configuration 

Level 

GRASP 

AS 

Index 

Configuration 

Level 

GRASP 

AS 

i 

2s 2 

1 So 

0.00000 

0.00000 

41 

2p3d 

3po 
r 2 

19.77978 

19.77212 

2 

2s2p 

3po 

*0 

1.41897 

1.41738 

42 

2p3d 

3po 

r l 

19.79099 

19.78275 

3 

2s2p 

3 po 
r l 

1.43387 

1.43437 

43 

2p3d 

3po 

r o 

19.79743 

19.78828 

4 

2s2p 

3po 

r 2 

1.46662 

1.46968 

44 

2p3p 

1 s 0 

19.94280 

19.92116 

5 

2s2p 

lpo 
r 1 

2.81633 

2.79949 

45 

2p3d 


20.02755 

20.01046 

6 

2p 2 

3 Po 

3.71907 

3.71638 

46 

2p3d 

1 po 
r l 

20.07597 

20.06076 

7 

2p 2 

3 Pi 

3.73652 

3.73517 

47 

2s4s 

3 Si 

22.54097 

22.53181 

8 

2p 2 

3 Po 

3.76518 

3.76832 

48 

2s4s 

1 s 0 

22.64271 

22.63203 

9 

2p 2 

3 d 2 

4.17504 

4.17207 

49 

2s4p 

3po 

r o 

22.79178 

22.78517 

10 

2p 2 

% 

5.17938 

5.16932 

50 

2s4p 

3 po 
r l 

22.79331 

22.78674 

11 

2s3s 

3 Si 

16.89542 

16.88436 

51 

2s4p 

3po 

r 2 

22.79723 

22.79057 

12 

2s3s 

% 

17.15693 

17.15202 

52 

2s4p 

lpo 

r l 

22.82350 

22.81370 

13 

2s3p 

lpo 
r 1 

17.52796 

17.51781 

53 

2s4d 

3 Di 

22.93591 

22.92974 

14 

2s3p 

3po 

*0 

17.55008 

17.54198 

54 

2s4d 

3 d 2 

22.93645 

22.93041 

15 

2s3p 

3 po 
r 1 

17.55555 

17.54746 

55 

2s4d 

3 d 3 

22.93741 

22.93144 

16 

2s3p 

3po 

r 2 

17.56323 

17.55522 

56 

2s4d 

3 d 2 

23.02777 

23.01752 

17 

2s3d 

3 Di 

17.90599 

17.89877 

57 

2s4f 

3 f° 

23.02937 

23.02096 

18 

2s3d 

3 d 2 

17.90756 

17.90073 

58 

2s4f 

3 f 3 

23.02961 

23.02133 

19 

2s3d 

3 d 3 

17.91022 

17.90367 

59 

2s4f 

3 po 
r 4 

23.03003 

23.02184 

20 

2s3d 

3 d 2 

18.17882 

18.16519 

60 

2s4f 

lpo 
r 3 

23.05392 

23.04587 

21 

2p3s 

3po 

*0 

18.69110 

18.68325 

61 

2p4s 

3po 

r o 

24.20663 

24.20036 

22 

2p3s 

3po 

18.70615 

18.69882 

62 

2p4s 

3 po 

24.21768 

24.21073 

23 

2p3s 

3po 
r 2 

18.74171 

18.73549 

63 

2p4s 

3po 

r 2 

24.25816 

24.25301 

24 

2p3s 

lpo 

r l 

19.00401 

18.99133 

64 

2p4s 

lpo 

r l 

24.32131 

24.30294 

25 

2p3p 

3 Pi 

19.08552 

19.07491 

65 

2p4p 

x Pi 

24.38599 

24.37956 

26 

2p3p 

3 Di 

19.16008 

19.14799 

66 

2p4p 

3 Di 

24.41706 

24.41031 

27 

2p3p 

3 d 2 

19.17515 

19.16319 

67 

2p4p 

3 d 2 

24.42038 

24.41402 

28 

2p3p 

3 d 3 

19.20844 

19.19688 

68 

2p4p 

3 d 3 

24.45399 

24.44816 

29 

2p3p 

3 Si 

19.32141 

19.30707 

69 

2p4p 

3 Si 

24.47906 

24.46575 

30 

2p3p 

3 Po 

19.38473 

19.37790 

70 

2p4p 

3 Po 

24.48248 

24.47470 

31 

2p3p 

3 Pi 

19.40025 

19.39402 

71 

2p4p 

3 Pi 

24.50993 

24.50041 

32 

2p3p 

3 p 2 

19.41709 

19.41265 

72 

2p4p 

3 P 2 

24.51423 

24.50850 

33 

2p3d 

3 f° 

19.46968 

19.46172 

73 

2p4d 

3 f° 

24.53225 

24.52677 

34 

2p3d 

3 f 3 ° 

19.49594 

19.48869 

74 

2p4d 

3 f 3 

24.55592 

24.55041 

35 

2p3d 

3 D 2 

19.50998 

19.50530 

75 

2p4d 

3 D 2 

24.56504 

24.56022 

36 

2p3d 

3 FS 

19.52225 

19.51643 

76 

2p4p 

3 d 2 

24.58233 

24.57285 

37 

2p3p 

3 D 2 

19.60698 

19.59904 

77 

2p4d 

3po 

r 4 

24.58565 

24.58108 

38 

2p3d 

3 d° 

19.68861 

19.68287 

78 

2p4d 

3 dj 

24.61082 

24.60428 

39 

2p3d 

3 d? 

19.69591 

19.69032 

79 

2p4d 

3 d° 

24.62033 

24.61401 

40 

2p3d 

3 D 3 ° 

19.70995 

19.70524 

80 

2p4f 

l f 3 

24.62916 

24.62052 


GRASP: |Aggarwal & Keenan|( 

2014a 

AS: |Fernandez-Menchero et al. 

(2014 


see for example levels of Ge XXIX in table 3 of|Aggarwal fe| 


Keenan (2014c I. 


Since two independent calculations for 5 Be-like ions 
with 13 ^ Z ^ 32 using the same R-matrix method (al¬ 
though in different approximations) and of similar complex¬ 
ity are now available, it is possible to make a detailed com¬ 
parison to assess their accuracy. This is important, given the 
large discrepancies recently noted for transitions of Fe XIV 


between our calculations with DARC (Aggarwal & Keenan 
|2014b| ) and those of |Liang et al.| ( |2010[ ) with ICFT. 


2 DETAILS OF CALCULATION 

In our work the fully relativistic GRASP (General-purpose 
Relativistic Atomic Structure Package) code is employed to 
determine the wave functions. There are several versions of 


this code, but all are based on the one originally developed 
by Grant et al. (19801, often referred to as GRASPO. The 
version used by us has been extensively revised by one of 
its authors (Dr. P. H. Norrington), is freely available at the 
website http://web.am.qub.ac.uk/DARC/ and yields com¬ 
parable results for energy levels and A-values as obtained 
with using other revisions. It is fully relativistic, based on 
the jj coupling scheme, and includes higher-order relativis¬ 
tic corrections arising from the Breit (magnetic) interaction 
and quantum electrodynamics effects (vacuum polarisation 
and Lamb shift). Additionally, the option of extended av¬ 
erage level has been adopted for all ions, under which a 
weighted (proportional to 2j+l) trace of the Hamiltonian 
matrix is minimised. 

For all ions the lowest 98 levels belonging to the 17 con¬ 
figurations (namely (Is 2 ) 2£2T / , 2131' and 2£A£') have been 
considered. As stated earlier, for the calculations of we 
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Table 2. Comparison of effective collision strengths (T) for resonance transitions of A1 X. a±6 = axlO ±b . 


Transition 

DARC (log T e , 

K) 

ICFT (log T e , 

K) 

i 

j 

4.30 

6.00 

7.30 

4.30 

6.00 

7.30 

i 

2 

1.751-2 

1.125-2 

2.333-3 

2.06-2 

1.09-2 

2.33-3 

i 

3 

5.841-2 

3.421-2 

7.510-3 

7.47-2 

3.31-2 

8.33-3 

i 

4 

1.127-1 

5.695-2 

1.173-2 

1.28-1 

5.43-2 

1.16-2 

i 

5 

1.128-0 

1.311-0 

1.940-0 

1.13-0 

1.30-0 

2.16-0 

i 

6 

3.925-4 

6.619-4 

1.485-4 

3.07-4 

5.40-4 

1.20-4 

i 

7 

1.021-3 

1.871-3 

4.029-4 

1.08-3 

1.59-3 

3.33-4 

i 

8 

1.975-3 

3.002-3 

6.714-4 

1.95-3 

2.69-3 

6.12-4 

i 

9 

1.448-2 

1.734-2 

1.286-2 

1.43-2 

1.69-2 

1.34-2 

i 

10 

5.105-3 

7.309-3 

3.964-3 

4.94-3 

7.47-3 

3.99-3 

i 

11 

6.103-2 

1.080-2 

1.291-3 

3.97-2 

9.33-3 

1.22-3 

i 

12 

9.534-2 

6.066-2 

6.852-2 

7.54-2 

5.98-2 

7.14-2 

i 

13 

4.994-2 

3.739-2 

1.327-1 

3.93-2 

3.66-2 

1.50-1 

i 

14 

5.753-3 

1.749-3 

2.783-4 

3.89-3 

1.58-3 

2.72-4 

i 

15 

1.965-2 

7.004-3 

8.022-3 

1.31-2 

5.89-3 

6.37-3 

i 

16 

3.012-2 

8.781-3 

1.393-3 

2.07-2 

7.92-3 

1.36-3 

i 

17 

1.232-2 

6.750-3 

1.304-3 

7.94-3 

6.13-3 

1.26-3 

i 

18 

2.039-2 

1.125-2 

2.179-3 

1.38-2 

1.02-2 

2.23-3 

i 

19 

2.843-2 

1.578-2 

3.047-3 

1.89-2 

1.43-2 

2.94-3 

i 

20 

5.683-2 

7.410-2 

1.424-1 

5.58-2 

7.11-2 

1.48-1 

i 

21 

8.029-4 

1.803-4 

1.379-5 

6.31-4 

1.53-4 

1.26-5 

i 

22 

2.980-3 

6.267-4 

1.037-4 

2.17-3 

4.80-4 

9.90-5 

i 

23 

4.853-3 

1.014-3 

7.485-5 

3.91-3 

7.54-4 

6.18-5 

i 

24 

1.127-2 

2.773-3 

4.427-3 

8.40-3 

2.56-3 

4.60-3 

i 

25 

3.740-3 

5.209-4 

8.951-5 

4.84-3 

4.82-4 

8.37-5 

i 

26 

2.525-3 

4.643-4 

6.364-5 

1.20-3 

3.85-4 

5.97-5 

i 

27 

4.160-3 

7.776-4 

1.041-4 

2.17-3 

6.40-4 

1.03-4 

i 

28 

4.976-3 

1.026-3 

1.347-4 

3.27-3 

8.73-4 

1.29-4 

i 

29 

5.235-3 

5.880-4 

5.978-5 

5.09-3 

6.35-4 

6.61-5 

i 

30 

7.920-4 

1.069-4 

1.160-5 

5.75-4 

6.49-5 

9.34-6 

i 

31 

4.201-3 

3.219-4 

2.905-5 

2.07-3 

1.96-4 

2.16-5 

i 

32 

6.179-3 

4.767-4 

5.373-5 

3.55-3 

3.12-4 

4.59-5 

i 

33 

1.855-3 

5.452-4 

8.664-5 

8.55-4 

4.88-4 

7.97-5 

i 

34 

1.673-3 

7.038-4 

9.523-5 

1.23-3 

6.45-4 

9.81-5 

i 

35 

1.645-3 

5.751-4 

1.294-4 

1.72-3 

5.51-4 

1.11-4 

i 

36 

1.892-3 

8.610-4 

1.157-4 

1.66-3 

8.18-4 

1.16-4 

i 

37 

1.967-3 

1.238-3 

1.164-3 

1.95-3 

1.23-3 

1.24-3 

i 

38 

1.888-4 

2.327-4 

8.316-5 

8.24-5 

1.48-4 

8.07-5 

i 

39 

3.312-4 

3.172-4 

3.247-5 

1.50-4 

2.21-4 

2.61-5 

i 

40 

4.338-4 

3.794-4 

3.830-5 

1.83-4 

2.57-4 

3.14-5 


adopted the DARC code. Specific details of the calculations 
are available in | Aggarwal fe Keenan| (|2014a[) for A1 X. 


garwal fc Keenan (2014c I for Cl XIV, K XVI and Ge XXIX, 


and Aggarwal & Keenan (2012b) for Ti XIX. Briefly, all 
partial waves with angular momentum J 40.5 have been 
considered, and to ensure convergence of II for all transitions 
and at all energies a “top-up” based on the Coulomb-Bcthe 
(B urgess fe Sheorey|[1974 l and geometric series approxima¬ 
tions was included for the allowed and forbidden transitions, 
respectively. Furthermore, values of SI were determined up 
to a wide range of energies up to 380 Ryd (A1 X), 660 Ryd 
(Cl XIV), 780 Ryd (K XVI), 1150 Ryd (Ti XIX) and 2500 
Ryd (Ge XXIX). 


For the subsequent calculations of T, resonances in a 
fine energy mesh (0.001 Ryd for most threshold regions) 
were resolved and averaged over a Maxwellian distribution 
of electron velocities. This distribution is commonly used 
and is appropriate for most astrophysical applications. The 
density and importance of resonances for Be-like ions can be 


appreciated from figs. 6-11 of Aggarwal & Keenan (2012b) 
for a few transitions of Ti XIX. Results for T were obtained 
over a wide range of electron temperatures (T e ) fully cover¬ 
ing that of their maximum fractional abundance in ionisa¬ 
tion equilibrium (Bryans, Landi & Savinj2009). Specifically, 
T were reported up to log T e = 7.2 (A1 X), 7.5 (Cl XIV), 
7.5 (K XVI), 7.7 (Ti XIX) and 7.8 K (Ge XXIX). 


Fernandez-Menchero et al. (2014) adopted the AS code 
of Badnell (199TJ) to calculate energy levels and A-values. 
Their calculations are comparatively larger as they included 
238 fine-structure levels of the (Is 2 ) 2(s,p) n£ (n = 3-7, l = 
0-4 for n ^5 and l = 0-2 for n = 6-7) configurations. For 
most Be-like ions considered, their semi-relativistic approach 
generally yields comparable results with GRASP for energy 
levels and A-values. For example, in Table 1 we list the two 
sets of energies for the lowest 80 levels of Al X, where dis¬ 
crepancies are less than 0.02 Ryd and level orderings are 
also the same. Since detailed comparisons of energy lev¬ 
els with available experimental and other theoretical results 
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Transition 

DARC (log T e , 

K) 

ICFT (log T e , 

K) 

i 

j 

4.30 

6.00 

7.30 

4.30 

6.00 

7.30 

i 

41 

1.054-3 

7.252-4 

1.069-4 

6.31-4 

7.78-4 

1.13-4 

i 

42 

6.700-4 

4.711-4 

8.345-5 

3.72-4 

4.84-4 

9.75-5 

i 

43 

2.164-4 

1.491-4 

2.203-5 

1.23-4 

1.61-4 

2.33-5 

i 

44 

1.531-3 

9.485-4 

5.759-4 

1.25-3 

1.11-3 

6.90-4 

i 

45 

9.098-4 

1.390-3 

1.151-3 

8.39-4 

1.44-3 

1.32-3 

i 

46 

3.728-3 

4.801-3 

9.630-3 

3.65-3 

4.87-3 

1.03-2 

i 

47 

9.284-3 

1.669-3 

3.098-4 

1.11-2 

1.94-3 

3.26-4 

i 

48 

1.830-2 

1.185-2 

1.384-2 

2.19-2 

1.15-2 

1.42-2 

i 

49 

1.342-3 

4.996-4 

9.613-5 

1.87-3 

5.41-4 

9.80-5 

i 

50 

4.530-3 

1.576-3 

5.461-4 

5.73-3 

1.70-3 

6.36-4 

i 

51 

6.416-3 

2.463-3 

4.768-4 

9.48-3 

2.71-3 

4.90-4 

i 

52 

1.088-2 

1.008-2 

3.006-2 

1.53-2 

9.68-3 

3.22-2 

i 

53 

3.009-3 

1.877-3 

4.081-4 

2.78-3 

1.84-3 

4.04-4 

i 

54 

4.942-3 

3.125-3 

6.814-4 

4.64-3 

3.07-3 

7.25-4 

i 

55 

6.894-3 

4.374-3 

9.513-4 

6.46-3 

4.30-3 

9.43-4 

i 

56 

1.425-2 

1.493-2 

2.732-2 

1.01-2 

1.30-2 

2.67-2 

i 

57 

1.685-3 

1.100-3 

1.973-4 

1.94-3 

1.27-3 

2.10-4 

i 

58 

2.613-3 

1.555-3 

2.778-4 

2.71-3 

1.78-3 

3.15-4 

i 

59 

3.032-3 

1.978-3 

3.547-4 

3.49-3 

2.29-3 

3.78-4 

i 

60 

6.258-3 

5.002-3 

7.684-3 

4.98-3 

5.34-3 

8.13-3 

i 

61 

7.022-5 

1.104-5 

1.956-6 

1.84-4 

2.84-5 

3.01-6 

i 

62 

2.462-4 

4.459-5 

1.894-5 

6.90-4 

1.03-4 

1.61-5 

i 

63 

3.067-4 

5.292-5 

9.511-6 

1.35-3 

1.42-4 

1.48-5 

i 

64 

5.487-4 

1.784-4 

1.698-4 

3.12-3 

2.94-4 

9.40-5 

i 

65 

1.626-4 

4.927-5 

1.409-5 

6.88-4 

1.16-4 

1.77-5 

i 

66 

1.310-4 

4.662-5 

1.255-5 

7.56-4 

1.20-4 

1.65-5 

i 

67 

2.124-4 

7.299-5 

1.818-5 

1.37-3 

2.26-4 

2.65-5 

i 

68 

2.643-4 

9.706-5 

2.190-5 

1.91-3 

3.08-4 

3.25-5 

i 

69 

1.368-4 

5.386-5 

1.198-5 

7.17-4 

1.43-4 

2.23-5 

i 

70 

2.675-5 

1.261-5 

5.294-6 

6.99-5 

3.10-5 

2.18-5 

i 

71 

1.219-4 

4.371-5 

9.004-6 

4.15-4 

7.54-5 

1.11-5 

i 

72 

1.125-4 

4.630-5 

1.164-5 

4.81-4 

7.88-5 

1.32-5 

i 

73 

1.545-4 

8.364-5 

2.041-5 

8.13-4 

1.52-4 

2.52-5 

i 

74 

1.877-4 

1.100-4 

2.396-5 

1.24-3 

2.06-4 

3.66-5 

i 

75 

1.657-4 

8.476-5 

2.622-5 

1.15-3 

1.51-4 

2.81-5 

i 

76 

2.350-4 

1.284-4 

1.143-4 

1.59-3 

2.43-4 

9.16-5 

i 

77 

2.456-4 

1.402-4 

2.784-5 

1.43-3 

2.56-4 

3.85-5 

i 

78 

7.624-5 

6.637-5 

8.803-5 

3.57-4 

9.66-5 

1.29-4 

i 

79 

1.111-4 

6.312-5 

1.416-5 

7.50-4 

1.11-4 

1.78-5 

i 

80 

4.524-5 

2.506-5 

8.451-6 

7.04-4 

8.66-5 

1.09-5 


DARC: Present results from the DARC code 
ICFT: Results of|Femandez-Menchero et al.||2014) 


have already been made by us and Fernandez-Menchero et| 
[ah| (2014]), we do not discuss these further. Similarly, we do 
not discuss the A-values, and rather focus on the more im¬ 
portant parameter, i.e. T. 

For the determination of fl, |Fernandez-Menchero et al.| 


(20141 obtained results for fine-structure transitions from 
the LS calculations ( |Berrington et aT7||1995[ ) through the 
ICFT approach. They considered a slightly larger range of 
partial waves (J ^ 44.5) than ourselves. However, they only 
included electron exchange up to J =11.5 and for the rest 
performed a no-exchange calculation for expediency. This 
approach sometimes leads to sudden changes (at the over¬ 
lap point) in the variation of If with J - see for example, 


table 6 of Aggarwal & Hibbert (1991). However, for most 


transitions it should not be a significant source of inaccu¬ 
racy. More importantly, Fernandez-Menchero et al. (2014) 


namely up to 90 Ryd (Al X), 165 Ryd (Cl XIV), 215 Ryd (K 
XVI), 300 Ryd (Ti XIX) and 680 Ryd (Ge XXIX), lower by 
a factor of ~ 4 compared to our work. Unfortunately, such 
energy ranges are insufficient ( Aggarwal &: Keenan|2008 1 to 
determine values of T at the high temperatures (up to ~ 
1.7xl0 9 K or equivalently ~ 10,600 Ryd) for which these 
authors reported results. They did include high energy con¬ 
tributions to fl from the suggested formulae of [Burgess fe| 
Tully (19921, but this approach, although computationally 


highly efficient, is perhaps a major source of inaccuracy, as 
discussed earlier by us for transitions of Fe XIV (Aggarwal] 
& Keenan 2014b I. Since we had already calculated values 


of S7 up to sufficiently high energies, there was no need for 
extrapolation to determine T for the ranges of T e reported 
by us. 


Fernandez-Menchero et al. (20141 resolved resonances 


performed their calculations of VL for limited energy ranges, 


in the threshold regions and averaged the values of over a 
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Figure 1 . Collision strengths for the 2s 2 1 Sq - 2p3p 3 Di (1-26) transition of A1 X. 


Maxwellian distribution of electron velocities to determine 
T, in a similar procedure to that employed by us. However, 
within thresholds they adopted a uniform mesh of 0.00001 
z 2 , where z is the reduced charge of the ion, i.e. Z-4. Con¬ 
sequently, with increasing Z their adopted energy mesh be¬ 
comes coarser. For example, for A1 X it is 0.00081 Ryd, but 
for Ge XXIX is 0.0078 Ryd. Furthermore, this energy mesh 
was adopted only for partial waves with J ^ 11.5 (i.e. for 
the exchange calculations only), while for higher J the mesh 
was coarser by a factor of 100, i.e. 0.001 z 2 . By contrast, we 
adopted a uniform mesh (~ 0.001 Ryd) for all partial waves 
and for all ions. 


3 EFFECTIVE COLLISION STRENGTHS 


Since Fernandez-Menchero et al. (20141 have not reported 
results for Q. no direct comparisons with our work are pos¬ 
sible. Therefore, we focus on a comparison of T, which are 
normally required for modelling applications. 


3.1 Resonance transitions 

In Table 2 we list both sets of results for T for the reso¬ 
nance transitions of Al X among the lowest 80 levels which 
have the same orderings in both calculations, and at three 
electron temperatures, i.e. log T e = 4.3, 6.0 and 7.3 K. The 
first and the third are the lowest and the highest common 


temperatures between the two calculations, whereas the sec¬ 
ond is the most relevant for modelling applications, because 
10 6 ' 1 K is the temperature at which Al X has its maximum 
abundance in ionisation equilibrium ( jBryans et al. [200 9). 

For most transitions listed in Table 2, the T of 
|Fernandez-Menchero et al. (2014| are larger, by up to a fac¬ 
tor of 15, at all temperatures, particularly the lower ones. 
However, for some transitions our results are higher by up 
to a factor of 2, such as 1-26/27/28 (i.e. 2s 2 1 So - 2p3p 
3 Di, 2 , 3 ). These transitions are forbidden in both the LS and 
jj coupling schemes, and hence resonances for these are sig¬ 
nificant (particularly towards the lower end of the energy 
range) as shown in Fig. 1 for 2s 2 1 So - 2p3p 3 Di (1-26). We 
note that even a slight shift in resonance positions can affect 
the calculated values of T, particularly at low(er) tempera¬ 
tures. However, as shown in Table 1 the two sets of energies 
obtained by the AS and GRASP codes are comparable for 
most levels, including 2p3p 3 Di, 2 , 3 . If these near threshold 
resonances are missing from the calculations of Fernandez-| 


Menchero et al. (2014) then their results for T will clearly 


be lower. 


Among the transitions for which Ticft are larger than 
Tdarc at log T e = 4.3 K, are 1-64 (2s 2 - 2p4s X P?) 

and 1-80 (2s 2 1 So - 2p4f 1 Fa), which are allowed and for¬ 
bidden, respectively. The 1-64 transition is weak with f = 
5xl0 -4 in our work (GRASP) and f = lxlO -4 from AS, 
and yet Ticft > Tdarc by a factor of ~ 6 at the lowest 
common temperature, while Ticft is smaller by a factor 
of 2 at log T e = 7.3 K. For this transition, resonances are 
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Figure 2. Ratio (R) of DARC and ICFT effective collision strengths for transitions among the lowest 80 levels of A1 X, at electron 
temperatures of (a) T e = 10 4 ' 3 K, (b) T e = 10 6 ° K and (c) T e = 10 7 ' 3 K. Negative R values indicate that Ticft > Tdarc- 


not prominent as expected. Similarly, for the 1-80 (forbid¬ 
den) transition we do not observe any significant resonances, 
although Fernandez-Menchero et al. (2014) may do so, be¬ 
cause their calculations include 238 levels in comparison to 
only 98 in ours. Therefore, we now focus our attention on the 
highest temperature of 10 7 ' 3 K (equivalent to ~ 126 Ryd) 
at which the contributions of resonances, if any, will not be 
appreciable. 

Among the transitions listed in Table 2, the one which 
shows the largest factor (4) by which Ticft is higher than 
Tdarc is 1-70 (2s 2 1 So - 2p4p 3 Po). In fact, the difference 
between the two calculations increases with increasing tem¬ 
perature (from a factor of 2.6 at T e = 10 4 ' 3 K to 4.1 at 10' 3 
K). However, this is a forbidden transition for which ft de¬ 
creases (or becomes nearly constant) with increasing energy 
- see table 4 of Aggarwal & Keenan (2014aI. Therefore, the 
behaviour of T in the calculations oflFernandcz-Mc ncherol 
|eFal. |2014| is difficult to understand for some transitions. 


3.2 All transitions 


The comparisons of T shown in Table 2 are for a very limited 
range of transitions. Therefore, in Fig. 2 (a, b and c) we 
compare the two sets of T for all 3160 transitions among 
the 80 levels of Al X at three temperatures of Table 2. It 
is clear from these figures that for a majority of transitions 
the T of Fernandez-Menchero et al. (2014) are significantly 
larger (by up to a factor of 30) at all temperatures. More 


specifically, for 69% and 42% of the transitions the two sets 
of T differ by over 20% at the lowest and the highest common 
temperatures of 10 4,3 and 10 7 ' 3 K, respectively. Also, for a 
majority of transitions the values of Ticft are larger. A 
similar comparison is found for transitions in other ions, 
namely Cl XIV, K XVI, Ti XIX and Ge XXIX. 


3.3 n = 2 transitions 

We now discuss the discrepancies between the two sets of T 
for specific range of transitions. Among the lowest 10 levels 
are 45 transitions belonging to the n = 2 configurations. At 
T e = 10 4 ' 3 K, for ~ 30% of transitions discrepancies are 
up to a factor of 2.5, largest for the 4-6: 2s2p 3 P 2 - 2p 2 
3 P 0 , which is forbidden. For half the transitions Tdarc > 
Ticft and for the other half Ticft > Tdarc, and there are 
no systematics because both types (allowed and forbidden) 
show the discrepancies. This may be due to the position of 
resonances, as discussed in section 3.1, because at T e = 10 6 '° 
and 10' 3 K both sets of T agree within ~20%. 


3.4 n = 3 transitions 

There are 630 transitions among the n = 3 configurations 
belonging to the levels between 11 and 46. At T e = 10 4 ' 3 
and 10 6 ' 0 K for about half the transitions (including both 
allowed and forbidden) Tdarc > Ticft, whereas a third 
of transitions have Ticft > Tdarc at T e = 10' 3 K. It is 
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difficult to explain these discrepancies at all temperatures, 
particularly for the allowed transitions. 


3.5 n = 4 transitions 

These transitions belong to levels 47 and higher and are 
561 in number. For a majority of transitions (about two 
third) Ticft > Tdarc at all temperatures. Again, these 
large discrepancies are difficult to understand, particularly 
at T e = 10 7 ' 3 K and for allowed transitions. Nevertheless, 
we will return to the source of discrepancies in section 4, 
apart from those already discussed in section 2. 


3.6 Strong allowed transitions 

We now discuss the strong allowed transitions, i.e. with f ^ 
0.1. Among the lowest 80 levels of A1 X given in Table 1, 
there are 79 such transitions. At the most relevant temper¬ 
ature of A1 X (i.e. 10® K) both sets of T agree within ~20% 
for all transitions, but for about a third Ticft > Tdarc 
at T e = 10 4 ' 3 and 10' 3 K. Discrepancies are comparatively 
larger at the lowest common temperature and many tran¬ 
sitions belong to the lowest 46 levels of the n = 2 and 3 
configurations. 


3.7 1—3 transition of P XII 


Finally, we note that the comparison shown by [Fernandez 


Menchero et al. (2014I in their fig. 5 with the T of 


Keenan 


(1988) for the 2s^ 1 So - 2s2p 3 Pj > (1-3) transition of P XII is 


incorrect. For this transition there is no appreciable differ¬ 


ence between the T interpolated by Keenan (19881 and those 
calculated by |Fernandez-Menchero et al. (2014). In fact, for 
this transition the interpolated values for other ions are un¬ 
derestimated (by up to a factor of two), rather than overesti¬ 
mated as shown by |Fernandez-Menchero et al. (20141. This 
can be seen from table 6 of Aggarwal & Keenan (2014a I for 


Al X, and tables 16 and 17 of Aggarwal & Keenan (2014c I for 


Cl XIV and K XVI. The error in calculating T has occurred 


by Fernandez-Menchero et al. (2014) because they have mis¬ 


takenly taken the coefficient al to be positive rather than 


negative (al = -0.026314), as given by Keenan (19881 in his 
table II. 


4 CONCLUSIONS 


In this paper we have compared two sets of electron impact 
excitation effective collision strength (T) for transitions in 
Be-like ions obtained by the R-matrix method in semi and 
fully relativistic approaches, i.e. ICFT and DARC. Both ap¬ 
proaches should provide comparable results for a majority 
of transitions. However, significant differences of up to more 
than an order of magnitude are noted for at least 50% of 
the transitions of all ions with 13 Z ^ 32, and over the 
entire range of electron temperature. In most cases the T 
from ICFT are significantly larger than those obtained with 
DARC. We believe the discrepancies have arisen mainly due 
to some compromises made by [Fernandez-Menchero et al.| 
(2014) in calculating and subsequently T, as noted in 
section 2. Similar large discrepancies between the two inde¬ 
pendent R-matrix approaches have also been noted in the 


past - see for example, Aggarwal & Keenan (2014b) for Fe 
XIV, Aggarwal & Keenan (2013b I and references therein 
for He-like and Aggarwal & Keenan (2012a I for Li-like ions. 
Therefore, it appears that the implementation of the ICFT 
approach, although computationally highly efficient (hence 
allowing data for many ions to be produced over relatively 
short periods), may not be completely robust. Indeed, this 
has also been confirmed by Storey, Sochi & Badnell (20141 in 
their calculations for O III, who noted that for some transi¬ 
tions the ICFT results can be significantly overestimated in 
comparison to the Breit-Pauli (or other similar approaches, 
such as DARC). Therefore, we recommend that the exci¬ 
tation rates reported by Fernandez-Menchero et al. (2014) 
should be used with caution and a re-examination of their 
results would be helpful. 

As already stated in section 1, assessing the accuracy of 
T is a difficult task ( Ag garwal fe Keenan|2013a l, mainly be¬ 
cause large calculations cannot be easily reproduced. How¬ 
ever, when two (or more) calculations of comparable com¬ 
plexity and with similar approaches (such as R-matrix) be¬ 
come available, the discrepancies observed are often striking, 
as noted above for He-like, Li-like and Be-like ions. There¬ 
fore, the true accuracy of any atomic data can only be as¬ 
sessed either by modelling applications or the comparisons, 
as shown in this work. However, since observational data are 
generally limited, a comparison with theoretical results is of¬ 
ten inconclusive. The large discrepancies in T observed here 
for a majority of transitions of many Be-like ions do create 
suspicion in the minds of users. Nevertheless, as discussed 
in our published papers, our accuracy assessment for a ma¬ 
jority of transitions in the temperature ranges of interest 
remains the same, i.e. ~20%. Although no accuracy assess¬ 
ment has been made by Fernandez-Menchero et al. (2014), 
their reported data would appear to be less reliable due to 
the compromises made in their calculations particularly in 
the chosen energy mesh and the energy ranges over which 
values of fl were calculated before extrapolation. However, 
we also stress that there is scope for improvement in our 
work (as there is for any calculation), especially for transi¬ 
tions involving levels of the n = 4 configurations. The inclu¬ 
sion of levels of the n ^ 5 configurations in the collisional 
calculations may improve the reported values of T due to 
the additional resonances arising. However, until then we 
believe that our reported results for radiative and excita¬ 
tion rates for transitions in Be-like ions are (probably) the 
most exhaustive and accurate available to date, and should 
be useful for the modelling of astrophysical plasmas. 
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